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0.  ABSTRACT 

A  finite  group  G  is  commonly  presented  by  a  set  of  elements  which 
generate  G.  We  argue y that  for  algorithmic  purposes  a  considerably  better 
presentation  for  a  fixed  group  G  is  given  by  random  generator  set  for  G: 

r 

a  set  of  random  elements  which  generate  G.  We  bound  .the  expected  number  of 
random  elements  required  to  generate  a  given  group  G. 

Our  main  results  are  probabilistic  algorithms  which  take  as  inputs  a 

•  :  »  •  \ 

— n.  ■  3 

random  generator  set  of  a  fixed  permutating  group  G  We  give  0  (n  log  n) 

expected  time  sequential  RAM  algorithms  for  testing  membership,  group 
inclusion  and  equality.  Our  bounds  hold  for  any  worse  case  input  groups? 
we  average  only  over  the  random  generators  representing  the  groups.  Our 
algorithms  are  two  orders  of  magnitude  faster  than  the  best  previous 
algorithms  for  these  group  theoretic  problems,  which  required  -fj(n^)  time 
even  if  given  random  generators. 

Furthermore,  we  show  that  in  the  case  the  input  group  is  a 
2-group  with  a  random  presentation,  than  those  group  theoretic  problems  can 
be  solved  by  a  parallel  RAM  in  O(log  n) ^  expected  time  using  n° ^ 
processors .  - 
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processors. 
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1.  INTRODUCTION 


I . I  Group  Inference  and  Representation  by  Random  Examples 

In  informal  mathematical  discourse,  we  often  make  use  of  examples  to 
illustrate  general  principles;  and  in  many  cases  this  suffices  to  convey 
the  essential  ideas.  For  instance,  Euclid  never  formally  stated  his  algorithm 
for  computing  GCD  but  instead  explained  it  by  a  series  of  example  computations. 
Cn  the  other  hand,  a  student  may  illustrate  comprehension  of  a  general 
principle,  say  Euclid's  algorithm,  by  producing  some  examples. 

The  f'.elds  of  Inductive  Inference  and  Combinatorial  Enumeration  concern 
the  dual  problems  of  inference  of  a  combinatorial  structure  from  examples, 
ar.d  generation  of  sample  examples  of  a  given  combinatorial  structure, 
respectively.  In  Section  2,  we  investigate  these  problems  when  the  combina¬ 
torial  structure  of  interest  is  a  finite  group,  and  the  samples  are  random, 
with  independent  uniform  distribution.  We  give  upper  bounds  on  the  number 
of  random  elements  of  G  required  to  generate  a  fixed  group  (generally,  the 
required  number  of  random  samples  is  a  logarithm  of  the  group's  order). 

As  an  interesting  example,  consider  the  group  of  all  permutations  of  the 
Rt'BIC's  cube.  Our  results  imply  this  group  can  be  generated  (with  high 
likelihood)  by  a  very  small  number  of  random  permutations.  Furthermore,  the 
results  allow  us  to  verify  (within  high  likelihood)  the  correctness  of  a 
"solution  method"  to  RUBIC's  cube  by  applying  the  "solution  method"  on  a  small 
number  random  example  permutations  of  RUBIC's  cube. 
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) . 2  Group  Theoretic  Problems 

The  fundamental  groups  problems  which  will  concern  us  are: 

(1)  Group  Membership  given  an  input  element  x,  and  a  group  G, 
test  x  £  G. 

(2)  Group  Inclusion:  given  groups  G,  H,  test  G  c  H. 

(3)  Group  Equality :  given  groups  G,  H  test  G=H. 

In  these  problems,  a  group  is  normally  assumed  to  be  finitely  presented. 
We  further  assume  that  a  group  is  randomly  presented  in  the  sense  defined 
below. 

Let  a  ixszdor.  generator  set  of  size  k  for  group  G  be  a  set  of  k 
random,  independently  chosen  elements  g  ,...,g  of  G,  with  the  condition 

-L  K 

that  q g.  generate  G.  Ke  say  <g,  , . . .  ,c,  >  is  a  rar.dcr.  presentation 
1  K  1  '  k 

c;'  G. 

A  probabilistic  alporith r  A  for  a  group  problem  takes  as  input  a 
random  generator  sets  cf  given  groups.  The  expected  time  complexity  of 

A  is  the  average  time  of  A  over  random  generator  sets,  given  worse  case 
groups . 

Section  2.~  describes  a  probabilistic  algorithm  for  constructing  a 
strong  generator  sequence  for  a  randomly  presented  group.  The  first  use  of 
such  probabilistic  constructions  (in  a  considerably  less  general  context 
was  due  to  iBabai,  79], 

We  also  discuss  in  Sections  2.3  and  2.6  known  algorithms  for  group 
membership  testing,  group  inclusion,  group  equality,  and  random  element 
generation  which  require  strong  generator  sequences. 


The  real  motivation  of  our  work,  and  our  strongest  results,  are  efficient 
probabilistic  algorithms  for  the  permutation  group  problems:  membership, 
inclusion,  and  equality. 

[Sim,  78]  first  gave  a  well  known  decision  algorithm  for  these  permutation 

group  problems  using  a  construction  known  as  the  Sim's  Table;  others  have  found 

it  to  be  very  efficient  in  practice.  However,  in  the  worse  case,  Sim's 

algorithms  were  exponential  time.  [Furst,  Hopcroft,  and  Luke,  81]  later 

0 

modified  Sim's  algorithm  to  yield  0(n  )  worse  case  time  bounds  for  these 
problems.  Still  further  work  by  [Jerrum,  82)  reduced  these  worse  case  time 
bounds  to  O(n^).  These  worse  case  time  bounds  seem  to  be  much  larger  than 
would  be  acceptable  in  practical  applications.  This  situation  motivated  us 
to  investigate  the  expected  time  complexity  of  permutation  group  problems. 

Our  main  results  are  C(n3  log  n)  expected  sequential  time  probabilistic 
algorithms  for  permutation  group  membership,  inclusion,  and  equality.  Our 
time  bounds  C(n3  log  n)  are  only  exceeded  with  probability  •  n~  '  for  some 
constant  a  1  that  can  be  set  arbitrarily  large.  In  comparison,  the  pre¬ 
viously  sighted  algorithms  for  these  permutation  group  problems  have  expected 
time  complexity  which  is  the  same  as  their  worse  case  complexitv  .  (n~). 

(Note :  Our  results  here  are  near  optimal  in  the  sense  that  we  can  show  that 

2 . 49 

further  decrease  in  our  sequential  time  bounds  below,  say  n  ,  would  imply 
improvement  in  the  sequential  time  bounds  of  best  known  algorithms  for  solving 
a  linear  system  of  size  n/2 * n/2  over  GF(2),  which  is  a  special  case  of 
permutation  2-group  membership  testing. ) 
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1 . k  Parallel  Algorithms  in  Group  Theory 

Ke  also  investigate  the  use  of  parallelism  for  group  theoretic  problems 
Our  goal  is  efficient  parallel  algorithms  requiring  polylog  time,  polynomial 
number  of  processors.  In  Section  4  we  give  efficient  parallel  algorithms 
for  thv  oriuits  and  block  systems  of  permutation  groups.  A  2-group  is  a  per¬ 
mutation  grout:  whose  elements  are  all  of  order  of  a  power  of  2;  they  arise 
naturally  since  a  2-crour  is  a  subgroup  of  automorphisms  of  binary  trees  and 
furthermore  the  automorphism  group  of  any  trivalent  graph  with  a  fixed 
vertex  is  a  2-group,  see  [Luks,  81].  We  give  in  Section  5  efficient  paral le¬ 
al  gorithnts  for  the  problems  of  membership,  inclusion  and  equality  for  2-grou 
Our  parallel  2-group-  membership  algorithm  makes  interesting  use  of  our  pro¬ 
babilistic  techniques :  it  takes  as  input  a  random  generator  set  cf  a  civ-;:. 
2-group,  and  constructs  from  them  a  tow-.r  :  f  C  v  1  c-  .*  r. '  subgre-uj  s  wit':.  h 
possible  to  do  efficient  parallel  membership  tests. 

The  parallel  complexity  of  the  permutation  group  membership  problem 
remains  open;  it  is  neither  known  to  be  log-space  complete,  for  deterministi 
polynomial  time,  nor  is  there  a  known  polylog  depth  algorithm,  for  group  me  mb 
ship.  Recently  [KcKer.zie  and  Cook,  83]  have  given  a  polylog  time  parallel 
algorithm  for  Abelian  permutation  group-  membership-.  Their  results,  and  our 
probabilistic  parallel  algorithm  for  2-group--  membership,  are  the  only 


positive  results  in  this  regard. 
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2.  PROBABILISTIC  ALGORITHMS  FOR  FINITE  GROUPS 


2. 1  Preliminary  Definitions 

Let  G  be  group  specified  by  [.resenting  a  finite  list  of  generators,  as  G 
<g  , . . .  ,g  > .  We  often  assume  G  has  a  finite  tower  of  subgroups 

X  -K 

G^  cG^c  ...  dG^  where  G^=G  and  G^  =  I  contains  only  the  identity  element. 

The  tower  has  height  h.  For  each  i=l,...,h  let  I ,  be  an  equivalence 

relation  such  that  Vx,y  £  G^  ,  x  =  .y  iff  y  ^x£g^.  The  blocks  of  each 

are  the  collection  of  cosets  of  G^  in  G^  ,  denoted  by  the  quotent 

3 .  , /G . .  Let  R.  be  a  complete  set  o'  oosez  representations  for  G.  ,/G.  , 
i-I  ii'  f  i-I  i 

i?~,  a  set  containin'?  exactly  one  element  from  each  coset  of  G.  _/G..  Since 

i-I  l 


=  (G,/G 


(G^/G^ )  . .  .  (G._  /G.  )  ,  the  sequence  of  sets  R^,...,R  is  called 


str.ng  so  :u^nee  of  generate rs  for  the  group  G,  with  respect  to  this  sub 


2 . 2  Elementary  Properties  of  RAND(G) 

We  will  let  RA::d(G)  denote  the  uniformly  distributed  random  variable 
giving  random  elements  of  G  (with  equal  probability).  If  x^,x^  are  rari^orn 
variables,  let  xg Jax2  ^  they  have  the  same  probability  distribution  function. 


LEMMA  2.1. 


*10  x  £  G,  Z'UW.  RAND  (G)  ^  x •  RAP D  ( o  )  . 


For  proof,  observe  that  the  function  f  (y)  =x*y  is  1-1. 


LEMMA  2.2. 


o"' coset  representatives  ~or  g/G ' .  ram D  ( g )  ®s RAND  ( R )  •RAND  (o’  ) 


P roof .  By  Lagrange's  Theorem,  each,  coset  of  G/G’  has  the  same  sice,  so 
RAND(G)  has  equal  probability  of  being  in  any  given  coset  of  G/G',  say  A. 

By  definition  RAND(R)  *  also  lias  this  property,  but  ;  roducting  an  element  of  A  by 
element  of  G'  keep  s  us  in  the  same  cosot  A.  Hence  RATIO (F)  1  RATIO (G'  )  has  the  same 


1 
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robafcility  of  being  in  A  as  does  RA!.TD(G).  By  definition,  RAh'D(G)  has  a 

liform  distribution  within  each  coset  A.  Let  x  be  the  unique  element  of 

A 


in  such  that  A  =  x 


Then  since- 


;  G ' | ,  we  must  have 


\ND(A)  ^  RAG'D  ( G 1  )  .  Hence  we  conclude  that  RAKD(R)  ^•RAIID(G')  is  also 

A 

niforr.lv  distributed  within  each  coset  of  G/G'  . 


, 3  Group  Inference  from  Random  Examples 

Let  G  be  a  fixed  finite  a  roup.  Let  L  =  ■ x  ,x0 , . . .  •  be  an  infinite 
ist  of  elements  chosen  independently  from  RAND (G) . 


HE  OR  EM  2.  1  .  0  <  •  1 ,  Frob  (G  =  <x,  ,  . .  .  ,x  >  )  ^  1 

1  m 


Z  +  (In  (log  C-)  'log  ( !/•:  ) ,  r:;e>v 


1  -  (1-0 


1/log  G 


roof  .  Let  J  =  i  €  •  1,2, . . .  •  x  .  €  <x  .  , . . .  ,x  .  ,  >  List  J  =  •  i  . 
-  1  1  1-1  O 


n  increasing  order.  With  probability  1,  this  gives  a  tower  of  subgroups 

~  ^  .  — '  •  .  •  ,  —  1  wnere  G  -  =  /  x  .  , .  •  •  ,  x  .  1*  ,  Bv  Laura  r,  u  c  1 : 

1  --Id;  .2  -s  1.  3  ' 

'  1  s 

<  G  /2~  and  J  <  loo  G  .  Hence- 


5  ~  J. 


3-1 


M"1  and  s:  I  rob  (i  >  j  +G)  <1/2  for  >  1 

'  s  s-1 


S  'S-1 


< 


s'io.Ol  ■  )  with  probabilil 


1  -  •  .  But  then  with 


rebabilitv  (1--  )"  >1--.  we-  have-  G=  <x  ,  ...,x  >  ii 

1  In 


I  » 

+  J2  (Js  -3S_!  -  V  <  ’-J'-  +  U 

s=2  s=2 


-  loq  (1/0 ! 
s  1 


J  ■  ♦  (In 


IDv 


)  log  (1 A 2  ) 


(In  (log  j  G  )  )  load/,,)  sir.c.  !j;  $  i0-  j  g 


a — an  a — x  w  i  i 
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POLYLOG  PARALLEL  TIME  ALGORITHMS  FOR  2-GROUPS 


A  finite  group  G  is  a  2-group  if  every  element  has  order  a  power  of 
Ke  will  show  that  any  2-qroup  has  a  certain  tower  of  h =  tlog  nj  subgroups 
i  — d  ->  .  ..-!  e>  =1.  Given  generators  for  each  subgroup  m  such  a 

:r,  we  can  test  membership  in  O(log  n)^  time  using  processors. 

:hermore,  if  G  is  given  by  random  presentation,  we  show  that  we  can 
struct  such  a  tower  in  0(loq  n)^  time  using  n0'^  processors. 


Parallel  Computation  of  a  Structure  Forest 


--crcu: 


can  be  decomposed  into  a  subgroup  of  natural  direct 


ducts  of  iterated  wreath  products.  Thus  there  is  a  structure  forcer  F 

err.;  lete  i  inary  trees  such  that  G  is  a  subgroup  of  the  natural  direct 

duct  cf  th>.  automorphism  groups  of  the  trees  of  F  .  In  particular,  each. 

t  cf  ’  is  called  a  err.-. mi .-.re  tv-.. •„  and  its  set  of  leaves  is  an  orbit  of 

Furthermore ,  let  FT  be  any  tree  or  subtree  of  F  which  is  not  a  leaf. 

G 

b  i-  tin  sc:  cf  leaves  of  ST  and  let  B^ ,  be  the  sets  of  leaves 

the  twe  imm.v  diute  subtrees  of  ST.  Then  we  require  that  -B  ,B?}  be  a 
lock  system  in  B. 

Suppose  we  are  given  generators  g  ,...,g  °f  2-aroup  G^S^.  By 
eating  first  the  G-orbit  algorithm,  and  then  executing  in  parallel  for 
a,b €  l,...,n:  the  G-block  algorithm  of  Section  4.2,  the  structure 

est  F  can  be  constructed  immediately  in  G(loq  n)  further  time  by 

G 

minaticn  of  each  G-block. 


MA  5.1.  eitK'.'jt  i v-.  ‘‘ov-.ct  F  o'  Z-orcuv  G  cs  can  be  cor.ctru 

’  G  n 

0(1) 

ref  a  a.'-.  1  t:r-  0(log  n)  ur :  kj  n  processors. 
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THEOREM  4.2.  .4  'i~'s  Table  jjk  be  constructed  fvc -m  a  rar.dcK  present. 


^  \  *  i.  .  c  C2  ^  4**  O  ^4  Y  «.  t  £_  4**  S  y  L-  }  *,  iT  vE'Z  tz.rc  0(n  log  n) 


COROLLARY  A.  I  . 


•rutatioyi  7rc:-.r2  in  s  .  liven  bu 

w  ‘  n 


"ft  -jtk  .* 


C(n  log  n) 


n 


Suppose  G  is  presented  as  <g  ,  For  distinct  a,  b€{l,...,n 

;t  us  construct  the  undirected  graph  with  vertex  set  il,....n}  ana  edqe- 
t  E  = { {a,b} }  U  { {g .  (a) ,g .  (b) } 1 1 ^ i  ^ k} . 

cL  f  D  1  X 

MMA  4.2.  [Atkinson,  75]  The  connected  component  of  ({l,...,n},E  ) 

a ,  b 

mzaining  a  is  the  smallest  G-block  containing  (a,b). 

Hence  finding  G-blocks  can  be  efficiently  reduced  to  undirected  graph 
>nnectivity,  which  can  apply  Lemma  4.1  to  get 

:MMA  4.3.  The  G-blcc'<c  can  be  computed  (in  the  worst  case)  in  time  oilog  n) 

3, 

n  k  processors. 

smark.  If  G  is  given  by  random  presentation,  G-blocks  can  be  found  in 
xpected  time  O(log  n)  using  n~  log  n  processors. 

•  3  Limited  Parallelism  for  General  Permutation  Group  Problems 

The  group  membership  algorithm  of  [Sims,  78]  was  improved  by  [Furst, 
operoft,  and  Luks,  80]  to  be  polynomial  time.  However,  it  appears  to  be 
nherently  a  sequential  algorithm;  and  can  not  apparently  be  speeded  up  to 
olylog  time  by  parallelization.  We  do  not  get  much  better  if  we  attempt  to 
irectly  parallelize  our  probabilistic  algorithms  for  general  permutation 
roups.  (But  our  low  processor  bounds  may  make  our  parallel  algorithms 
ore  practical.)  We  first  observe: 

EMMA  4.4.  Given  a  Sims  Table  for  a  permutation  croup  in  s  .  we  can  execute 

- r  n 

membership  test:  x  €  G?  in  time  o(n)  using  n  i  re  cess,  rs. 

Applying  Theorem  4.1  to  our  algorithm  of  Section  3.3  for  Sim's 


able  construction,  we  get: 


it.  PARALLEL  ALGORITHMS  FOR  PERMUTATION  GROUPS 

it.  1  The  Parallel  Machine  Model:  Known  Results 

We  will  assume  the  concurrent  read,  concurrent  write  parallel  RAM 
described  in  [Shiloach  and  Viskin,  82]. 

LEMMA  4.1  [Shiloach  and  Viskin,  82]  and  [Viskin  and  Tarjan,  84].  Given  an 
undirected  graph  of  n  vertices  arid  m  edges,  the  connected  components, 
a  spanning  forest,  and  a  preorder  of  each  tree  in  the  forest  can  all  be 
computed  in  time  o(log  n)  and  n+m  processors, 

4 . 2  Parallel  Computation  of  Orbits  and  Blocks  of  Permutation  Groups 

Let  G  ^ Sn  be  permutation  group  over  fl,...,n}.  It  follows  immediately 
from  Proposition  3.1,  Corollary  3.1,  and  Lemma  4.1  that 

THEOREM  4.1.  We  can  compute  the  orbits  of  G=  <g  ,...,g  >  in  time  O(log  n) 
using  in  the  worst  case  0(nk)  processors.  Furthermore,  if  G  is  given  by 
random  presentation,  we  can  compute  the  G -orbits  in  time  o(log  n)  using 
n  log  n  processors,  with  likelihood  1-n  for  any  sufficiently  large 
:nt  Ot  >  1  - 

G  acts  transitively  on  its  orbits.  G  is  transitive  if  G  has  only 
one  orbit.  Suppose  G  is  transitive.  A  nonempty  set  BC£{l,...,n}  is  a 
G-block  if  Vr,r'  EH,  r(B)=TT'(B)  or  rr  (B)  fl  ti  '(B)  =  0.  If  so,  then 
i ^  (B)  | -  €  G>  is  a  G -block  system,  and  group  G  acts  transitively  on  each 
of  the  blocks  of  the  system.  If  there  are  no  G-blocks  of  size  at  least  two, 
then  G  acts  primitively .  A  G-block  system  is  minimal  if  G  acts 


primitively  on  each  block. 


COROLLARY  3-2. 


For  (verse  case)  pcrr.uiaiicr.  groups  in  given  by  rancor 

presentation ,  permutation  group  rerberskip ,  inclusion  and  equality  can  all 
be  done  in  expected  time  0(n"^  log  n) .  Forth errr.o re,  these  bounds  hold  with 
probability  1-n  a  for  any  sufficient  '.y  large  cons- mi  a  ->  1. 


We  can  easily  compute  a  spanning  tree  T\  and  its  preorder  in  sequential 

time  0(n  log  n)  by  depth  first  search.  It  is  easily  to  verify  that  for 

each  j  £  V.  ,  r.  .  is  an  element  of  G.  ,  such  that  r.  .  (i)=i.  Hence 
i  1,3  i-l  i,D 

R.  =  {r.  . | j £  V . }  is  a  complete  set  of  coset  representatives  for  G.  /G ., 

l  i  /  3  D  -I-  i 

as  required.  Furthermore  L.  =  (r.  ^  .  •tt|tt£l.  -  F.  ,}  is  a  list  of 

i  i,ir(i)  i-l  i-l 

(n-i-1)  (ot+1)  c  log  n  elements  of  RAND(G.  Thus  F  . .  . .  ,  R  are  a  Sim's 

i-l  1  n 

Table  with  probability  ^1  -  n  a. 

The  most  costly  step  of  each  iteration  is  [10],  which  takes  time 

2  3 

C(n  log  n) .  Since  there  are  n  iterations,  the  total  time  is  0(n  log  n)  .o 


3.4  Sol ution  of  Permutation  Group  Problems  Utilizing  the  Sim's  Table 


The  three  lemmas  below  follow  immediately  from  the  discussion  in 
Section  2.  ’•>. 


LEMMA  3-2.  Giver,  the  Sir’s  Table  o'  a  reimutaticn  r.rcv.o  G  e  s  ,  tend  an 

•  -  1  —  n 

? 

input  x€  s  ,  Sir’s  "crbcrsktr  test:  x€  G?  takes  seauential  tore  c(n“). 
n 

LEMMA  3-3-  Given  permutation  group  G  =  <g  ,...,g  >  and  the  Sir's 

Sable  for  permutation  croup  g_  where  G, ,  G_  c  s  ,  then  the  a roue 
inclusion  test:  G1  c  g2?  takes  sequential  tire  Olnk^). 

Proof.  It  suffices  to  test  g^  €  G2  for  i=l,...,k  . 

Lemma  3.4.  Given  permutation  groups  G^  =  <g^,...,gk  >  and 

G  =  <h  ,...,h  >  in  s  and  their  Sim's  Tables ,  then  we  can  test  group 

z  x.  k  _  n 


'quality:  G  =  G1?  in  sequential  time  0(n(k^+k2) 


As  an  immediate  consequence  of  the  above  Lemmas  3. 2-3. 4,  Theorem  2.1, 


and  Corollary  3.1,  we  have: 
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begin 

Let  Lq  be  a  list  of  m’  random  elements  independently  drawn  from  RANT'(G) 


for  i 

=  1 , .  .  .  ,  n  do 

begin 

[1] 

Let  F,  be  the  set  consisting  of  the  first  (a+l)c 

log 

n 

elements  of  L.  , 
l-l 

12} 

Compute  the  connected  comp>onents  of  the  graph 

<:i . n}'  W  - 

1-1 

E  )  . 

TT 

[3] 

Let  V,  be  the  connected  component  containing  i. 

14} 

Compute  a  spanning  tree 

of  component  \h  rooted 

at 

i  . 

15] 

Label  each  edge  e  £ 

with  a  permutation  £(e)  =  ~ 

where  "6  F. 

1 

and  such  that  e£  E_. 

[6] 

Let  r.  .  be  the  identity  permutation 

17} 

In  a  preorder  traversal 

of  tree  Th ,  compute  for  each 

j  £  Vi  -  {i  } 

the  permutation  r 

.  =  r.  .  ,  •  H  (i  ’  ,  j )  where  j  ’ 
,3  i.3  ' 

is 

the 

parent  of  j . 

[S] 

K.  *-  {r.  .  |  j  £  V.  }. 

i  i,D  i 

[9] 

L  +■  0. 

[10] 

for  each  n £  L.  -  F.  , 

-  l-l  l-l 

do 

add  r  ^  .  *TT  to 

l  ,V  (l) 

L.  . 

l 

end 

return  the  Sim’s  Tabic  R,  ,  .  .  .  ,  R  . 

-  1  n 

For  proof  of  this  algorithm,  let  us  assume  inductively  for  some  i  >  1 

that  Lj_  j  is  a  of  (n-i)  (a+l)c  log  n  elements  of  RAND(G^_^) •  By 

"Ot 

Theorem  3.1,  \L  is  the  G^-orbit  containing  i,  with  probability  P  1-  n 
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Let  G*  =  ( -  -  -  (G  •*-)Tr  •  X  (G*)  is  just  the  number  of  connected 

12  ac  log  n 

components  of  graph  ({l,...,n),  li  ...  uE  )}  which  are  proper 

1  ac  log  n 

subsets  of  orbits  of  G. 

Lemma  3.1  implies  for  any  sufficiently  large  constant  a  above  c, 

that  Prob (I (G*)  =  0)  ^  1  -  -qj  ,  which  proves  our  Theorem  3.1.  O 

n 

3- 3  Constructing  a  Sim's  Table  in  0(n^  log  n)  Expected  Sequential  Time 

Fix  a  permutation  group  G  c  Sn  over  points  {l, . . . ,n}.  For  i =  1, . . . ,n, 

let  G^  be  the  subgroup  of  G  fixing  points  l,...,i.  The  resulting  tower 

G  =  G  ^  G  o  .  , .  a  G  =1  has  height  h  =  n  and  is  called  the  point  stabilizing 

tower  of  G-  A  strong  sequence  of  generators  R^ , .  .  . . R^  for  this  point 

stabilizing  tower  is  called  a  Sim's  Table. 

Unfortunately,  it  is  very  expensive  to  construct  the  Sim's  Table  by 

known  techniques.  [Furst ,  Hopcroft  and  Luks,  80)  give  the  first  polynomial 

0 

time  algorithm,  running  in  sequential  time  0(n  ) .  [Jerrum,  82]  improved 
this  time  to  the  best  known  worse-case  bound  of  O(n^). 

THEOREM  3-2.  Given  a  random  presentation  of  a  given  ( worse  ease)  permutation 

croup  G  c  s  ,  we  can  construct  a  Sim's  Table  in  expected  sequential  time 
w  —  n 

3  ■“  CX 

O (n  log  n).  Further,  these  bounds  hold  with  probability  >1  -  n  ,  for  a  su.ff 

ciently  large  constant  a  > 1  which  can  be  set  arbitrary ly  large. 

P roof .  We  will  fix  m'  =  (a+l)cn  log  n,  where  c  >1  is  a  sufficiently  large 


constant  and  a^c  is  an  arbitrarily  large  constant. 


algorithm  of  [Hopcroft  and  Tarjan,73]. 


COROLLARY  3.1.  We  can  compute  the  G-crbits  from  G(logn)  element c  of  R/UL(G)  in 
0(n  log  n)  sequential  time,  with  error  probability  <n  a  for  ar.y  suffi¬ 
ciently  larae  constant  a>l. 


Now  we 

prove 

Theorem  3. 

1 . 

Fix  a  permutation 

TT,  £  G.  Let  B  .  ,  .  .  .  ,B_ 

i  v1 

be  the  orbits  of 

the  group- 

<7 

1 

> .  For  each  tt  £  G , 

let  Cp_  (7)  be  the  per- 
"l 

mutation  ~ 

>  e  s 

k 

such  that 

Vj 

=  1 - k,  7-  (j)  = 

j  ’  iff  3i  s.t.  7  (i) £  E 

and  i  £  B_ 

Thus  ip.,  (t 

■) 

is  derived  from  7 

by  collapsing  each  orbit 

*1 

of  tt  to  a  single  point.  Let  G_  =  {tp_  (r)  !tt£  g).  Let  T  (G)  be  the  size 
1  '1  "1 


the-  set  ; i  (i  )  *  i  for  some  -  £  C-  • . 


1 


,k 

1* 

of 


LEMMA  3-1.  I -  77,  £  RAND ( G)  ,  then  prob(T(G  )  <T(G)/2)  . 

P roof .  Suppose  not,  then 


i  {  (77 ,  i  ) ;  r  £  G,  r(i)  jt  i  }  }  <  j  |g|t(G) 


By  the  pidgeon  hole  principle,  3iQ  such  that  HqUq)  f  iQ  for  some 
£  G  but 

I  {77  e  G  I  77  (iQ)  =  iQ  ) !  >  ^ 

But  this  implies  that  the  proper  subgroup  {7  €  G ' T  (iQ)  =  iQ-  of  G  has  order 
greater  than  |g|/2,  a  contradiction  with  Lagrange's  Theorem. 
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3.  PROBABILISTIC  ALGORITHMS  FOR  PERMUTATION  GROUPS 

Our  main  results  concern  permutation  group  problems.  Let  denote 

the  group  of  all  permutations  over  n  points. 


3 . 1  Inference  of  a  Permutation  Group  from  Random  Examples 

Theorem  2.1  implies  that  for  a  fixed  permutation  group  G  c  s  , 

2  2 

log  (n!  )  +O(log  n)  ^n  log  n-n  log  e+0(log  n)  independently  chosen 
permutations  from  RAND (G)  suffice  to  generate  G  with  probability  at 
least  1-n  U  for  any  sufficiently  large  constant  a  > 1. 


3 . 2  Computing  Orbits  in  Expected  Sequential  Time  0(n  log  n) 

Let  G  c  be  a  permutation  group  over  {l,...,n}.  The  G -orbit  of 
i€{l,...,n}  is  {tt  (i  )  ;  m  E  g)  .  Note  that  the  G-orbits  partition  {l,...,n;.  Le 
E_  =  {  (i  ,  j  )  '  tt  (i )  =  j  or  T(j)  =  i}. 

It  is  obvious  that  if  G  =  <g  , , . . ,g  >  then 

PROPOSITION  3-1-  The  G-orbits  are  the  connected  components  of 

({l, .: . , n } ,  U  .  E  )  . 

i=l , . . . ,k  g . 

l 

We  will  show: 


THEOREM  3.1.  For  all  a  above  a  constant  c>0, 
tt  €  RAND(G),  *  he k  v;  th  probability  at  least 

connected  components  of  the  graph  ({l,...,n},  E 


if  d  =  ac  log  n 


n 


TT  U  •  •  ■  °  Er  )  • 
1  d 


and  if 


We  can  compute  the  connected  components  of  a  graph  of  n  vertices  and 


cn  log  n  edges  in  sequential  time  Otn  log  n)  using  the  depth  first  search 
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2 . 6  Random  Element  Generation  from  a  Strong  Sequence  of  Generators 
Again  suppose  we  have  a  strong  sequence  of  generators 
for  group  G  with  respect  to  its  subgroup  tower  G  =  G^o  . .  .=>  G^  =  I. 
We  can  compute  RAND(G)  by  a  simple  algorithm  described  in  [Hoffman,  82] 


begin 


for  i  =  l,...,h  let  x.  be  a  random  element  of  R, 

-  l  l 

return  x, . . . x 

-  1  rc 


To  justify  this  algorithm,  observe  that  Lemma  2.2  implies: 

LEMMA  2.5.  RAND (G)  «  RAND(R;  )_1  . . .  -RAND (Ri ) _1 . 

Remark.  It  is  interesting  to  observe  that  a  random  element  of  G  can  be 
generated  by  this  method  in  parallel  by  a  binaiy  product  tree  of  depth  O(log  h) 


0 


I 


r~ 

> 


i 


i 


i 


Hence  we  have  total  failure  probability  at  most  hw(l-l/w)m  <  c .  If  there 


is  no  failure,  then  each  R^  is  a  complete  set  of  coset  representatives  for 


^/G.  .  Hence  are  the  strong  generators  for  this  subgroup 

tower.  □ 


2 . 3  Group  Membership  Inclusion,  and  Equality  from  Strong  Generators 

Let  G  =  G  r>G  . . .  oG,  =  I  be  a  subgroup  tower.  Let  R,  , . . .  ,R,  be  a 
0  i  h  1  h 

strong  sequence  of  generators  of  G  computed  in  the  previous  subsections. 

We  also  assume  that  for  any  i  and  x  €  G^  ^  we  can  effectively  find  the 
coset  representive  y€R.  such  that  x  =.  y  (ie,  so  y  ^x£G.). 

We  now  describe  Sim's  algorithm  for  membership,  in  the  general  context 
of  finite  groups. 

Given  an  input  x  and  such  a  strong  sequence  of  generators,  the  Sim's 
membership  algorithm  is: 
for  i  = 1  to  h  do 

if  3y  €  R.  s.t.  x  y  then  x  •>- y 
else  return  ("x  is  not  a  member  of  G" ) 

return  ("x  is  a  member  of  G")  . 

Remark.  Sim's  membership  algorithm  seems  inherently  sequential,  since  the 
parallel  time  for  its  execution  is  at  least  f.  (h),  where  h  is  the  height 
of  the  subgroup  tower. 

Also  observe  that  given  another  finitely  presented  group  G'  =  <g  ,...,g^, 
and  a  strong  generator  sequence  for  G,  we  can  test  group  inclusion  G'^G  by 
by  simply  testing  g^  £ G  for  i=l,...,k’-  Furthermore,  given  strong  generator 
sequences  for  finitely  presented  groups  G=  <h^,...,h^>  and 
G'  =  <g^, . . . ,gk,  > ,  we  can  test  G  =  G'  by  varifying  h^  €  G'  for  each 

k' . 


i  = 1, . . . ,k  and  g .  €  G  for  each  j  =  1, . . . , 


Let  Lq  be  a  list  of  m  elements  independently  drawn  from  RAND (C j  . 


for  i  = 1, . . . ,h  do 
begin 


R.  *  i; 

1 


for  each  A  €  G.  _/G.  do 
-  l-l  l  — 


begin 


Let  L .  .  be  the  list  of  elements  of  L.  ,  in  A 
l , A  l-l 


if  L.  „  is  not  ennty  than 
—  i,A 


begin 


choose  and  delete  a  random  element  r  fr 

A 


add  r  to  R. 

A  l 


for  each  remaining  element  xtL.  .  do 

-  i,/\ 


-1 


add  r  x  to  L. 


end 


end 


end 


return  strong  sequence  of  generators  P.^  ,  .  .  .  ,  R^ 


end 


THEOREM  2.2.  If  m>hw+  (log  t)/log  (hw(l-l/w))  then  the  ale  or;,  ter:  out; r;<r 
a  et  K'  e- j  -  •  ••  ;*  generators  of  G  with  error  probability  <  e  . 

Proof.  We  inductively  assume  that  on  the  i-th  iteration,  we  have  had  no 

failure  and  L .  ,  is  a  set  of  at  least  m  -  (i  -l)w  elements  independently 
i-l 

chosen  from  RAND (G .  Then  Lemma  2.4  implies  each  element  of  L.  is 

i-l  a 

independently  distributed  as  RAND(G^),  and  clearly  |  |  ^=|l^_^|  +  w^m-iw. 

By  Lemma  2.3,  the  probability  of  failure  at  stage  i  is  at  most  w(l-l/w) 
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2.4  Constructing  Strong  Generators  of  a  Group  from  Random  Examples 


Let  G  be  a  finite  group  with  subgroup  G'cG.  We  now  wish  to  discover 
a  complete  set  of  coset  representatives  of  G/G1.  Again  let  L  be  a  random 
list  of  m=  [  L  j  elements  chosen  from  independent  craws  from  RAND  (G)  . 

Let  E  be  the  event  that  L  contains  at  least  one  representative  of  each 
coset  in  G/G' . 


LEMMA  2.3.  If  G/G'  has  w  aoseis,  then  Prob (E)  >  1  -  w (1  -  1/w)  . 


Proof .  Consider  a  given  coset,  say  A,  of  G/G'.  Since  a  random  element 
of  G  has  equal  likelihood  to  be  any  coset  of  G/G', 

Frob  (LAA  =  : )  =  (1  -  l/w)m.  Hence  Prob  (not  E)^w  Prob  (LAA  =  :)  ^  w  ( 1  -  1/w) n. 


Let  us  assume  event  E.  For  each  coset  a£g/G',  fix  r  to  be  a  random 

A 

element  of  LaA,  and  let  R  =  {r  ;A€G/G'}.  For  each  x£G,  let  f  (x)=r  ^x 

A  R  A 

where  x  £  A. 


LEMMA  2.4.  f  (RAND (G) )  »  RAND(G'). 

R 

Proof.  Clearly  f„(x)  €G'  for  each  x£G.  By  definition,  f_(RAND(G))  has 
1  ~  R  R 

a  distribution  function  identical  to  r  ^ (RAND (A) )  for  randoml%-  chosen 

A 

A  £  G/G ' .  Also  since  r^A=G'  for  each  A£G/G',  we  have 

r*"1  (RAND  (A)  )  RiRAND  (G' )  .  Hence  f  (RAND (G) )  « r” 1  RAND  (A)  « RAND  (G '  )  . 

A  R  A 

We  now  assume  group  G  has  subgroup  tower  G  =  G^::;G^::3  . . .  G^  =  I  of 

height  h  and  width  w  = max  |  G^  /G^ j .  The  following  procedure  constructs  a 

i 

list  of  strong  generators  for  G  with  respect  to  this  subgroup  tower. 
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Tk 


The  structure  forest  F  aids  us  in  designing  an  efficient  parallel 

G 

algorithm  for  2-group  membership  testing. 

5.2  Root  Actions  are  Linear 

We  wish  to  reduce  2-group  membership  testing  to  solving  linear  ver 

GF(2)[x].  We  will  fix  throughout  Sections  5.2  and  5.3  the  2- 

Let  a  ,...,a  be  the  roots  of  structure  forest  F  .  Let  as t 

J.  2T  G 

on  root  a^  if  a^  is  not  a  leaf  and  it  permutates  the  two  children  of  a. 

and  otherwise  r  stabilizes  a^  For  any  -  £ Sn,  let  A  (-)  =  (A  (-),... ,A  (-) )T 

where  A.  (7)  =1  if  -  acts  on  the  root  a,  and  otherwise  let  A.  (-)  =0. 

3  3  3 

It  is  easy  to  verify: 

PROPOSITION  5.1.  V-  7  £s  ,  A(tt  •-  )  =  A(7  •-  )  . 

1  2  n  12  2  1 

Proof .  A(Vt2)  =A(7i)  +A(72)  =a(72)  +A(1)  -a<V*  ).  a 

Thus  the  permutations  act  commutitively  on  the  roots  of  the  structure 
forest . 

To  avoid  repeated  use  of  the  transpose  symbol,  we  will  simply  let 
k 

x  £  {0,1;  denote  a  column  vector  of  booleans.  Let  M  be  the  r  x  k  boolean  matrix 

such  that  Vi.l^i^k,  the  i-th  column  of  M  is  A(g  ).  Let  A(G)  =-A(~)  -  €  G ; . 

LEMMA  5.2.  A (G)  is  the  linear  space  (Mxjx  £ {0,1}  )  over  GF(2)[xJ. 

Proof .  Suppose  r,  £g.  Let  ”  =  g .  ...  g.  be  a  factorization  of  "  to  a 

]1  3s 


• 

product  of  its 

generators. 

For 

i  —  1 , . . . , k, 

let  x . 

1 

be  the  residue  mod  two 

of  |{t|jt=i} 

I  (ie. 

x.  =0 
1 

if 

g ,  occurs 

1 

an  even  number  of  times  in 

this 

factorization 

X 

of  7  , 

X 

and  x 

.  =  1 
1 

otherwise ) . 

Then  by 

Proposition  5.1, 

* 

A (~ )  =A(gxX  .. 

k, 

*  gk  > 

But 

the 

j-th  element 

of  AC") 

is  the  mod  two 

v 

sum  Z  x. 

1 

where 

the  sum 

is 

taken  over  just  those 

i  such  that  g. 

1 

acts 

;  ' 

on  a • .  Hence 

by  definition 

,  Mx  =  A  (7  )  . 
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Let  G^3  be  the  subgroup  of  G  consisting  of  the  permutations  that  fix 

(1 )  ,  t 

the  roots  a  , . . . ,a  ,  ie,  G  =  {tt  €  G  |  A  (ti  )  =  (0, . . .  ,0)  }. 

X  X 

LEMMA  5.3*  it  €  G  iff  3x£{0,l}k  s.t.  Mx=Ad)  and  (g  1  ...  g  k)  ^h€g^\ 

xi  xv  (1)  xi 

Proof.  If  Mx  =  A  ( ~ )  and  (q„  . . .  q,  )  it  £  G '  then  since  q.  . .  .  q,  £G, 

-  1  k  1  k 

we  have  tt  £G.  Suppose  on  the  other  hand  that  r  £ G.  By  Lemma  5.2, 

k  x  x 

3x  £  { 0, 1;  s.t.  Mx  =  A  H  )  and  so  A  (g^  ...  g^  )  =  A  (z ) ,  hence 

X1  Xk  -1  T  X1  Xk  -1  (1) 

A((91  *’•  gk  '  =  (° . 0)  50  (91  gk  }  r€G  •  D 

5 . 3  2-Group  Membership  Testing  Given  a  Block  Structure  Tower 

Let  the  block-s tructu re  tower  of  G  be  the  sequence  of  subgroups 

G  =  G^°^  3...  3G^'  =1  where  G^  is  the  subgroup  of  G  containing 

only  permutations  that  fix  all  the  nodes  of  depth  <  i  in  the  structure 

forest  F  .  Since  the  depth  of  F_  is  at  most  tlog  nj ,  this  tower  has 
G  G 

depth  n  ^  l log  nj  . 


THEOREM  5.1.  Suppose  we  have  generators  for  the  clock  structure  tower 
G  =  G(0)  =>G(1)  =>...  =>G(h)  -  I  of  the  2-group  GEsn-  2 hen  we  can  test 
in  G  in  time  0(log  n)"3  using  n° ^  processors. 


Proof.  Suppose  we  are  input  some  permutation  r  £  s  .  We  first  must  determine 

-  n 

the  existence  of  a  solution  x£{0,l}  of  the  linear  equation  Kx=A(~) 


defined  in  Lemma  5.2.  If  no  solution  exists,  we  reject 


this  solution  x  to  reduce  the  problem  to  membership  testing  in  G 


Otherwise  we  use 
(1) 


In 


particular,  by  Lemma  5.3,  we  need  only  test  if  (g  ...  g  )  z£g  This 


V-1 

3k 


(1) 


is  done  by  recursive  application  of  the  membership  test.  Since  h^  llog  nj , 

at  most  log  n  stages  suffice.  Each  stage  thus  requires  solution  of  a  linear 

system  over  GF(2)[x]  of  size  at  most  n*n,  which  can  be  done  in  time  O(log  n)‘ 
0(1) 


using  n 


by  the  parallel  algorithm  of  [Borodin,  van  zer  Gothen,  and  Hopcrc-ft, 
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82] 


Thus  the  total  time  is  0(log  n)  using  n 


0(1) 


processors. 


2 

Remark.  These  time  bounds  can  be  decreased  to  0(log  n)  if  the  required 
matrix  inverses  are  precomputed,  (so  that  each  stage  requires  only  matrix 
multiplications,  which  takes  only  0(log  n)  time  using  n^ ^  processors). 


5. k  Poly  log  Time  Construction  of  the  Block-Structure  Tower  from  a  Random 
Presentat i on 

Let  L r  }  be  a  list  of  m  permutations  independently  chosen 
1  m  o 

o 

from  RAND(G).  Fix  m  = m^/log  n.  Let  Y  be  the  linear  space  over  GF(2)[x] 

aenerated  by  A  (",),...,  A  (r  ).  Since  A(G)  is  a  qrouo  of  size  at  most  2r^2 
1  m  *  , 

by  Theorem  2.1  we  have 


LEMMA  5.^4.  Prob  (Y  =  A (G)  )  ^  1  -  e  if  m  >  r  -  (In  (r)  )  log  (1/c^)  ,  inhere 

Note  that  since  the  total  number  of  nodes  of  the  structure  fore 

2n,  we  can  bound  the  probability  of  error  to  be  at  most  £  using 

r  =  2n  +  (In (2n) ) log (1/s^)  random  elements  of  G.  For  example,  if 

■_  =  n  ,  then  it  suffices  to  let  m  =  2n+o(n). 

Let  M'  be  the  r  ■  m  boolean  matrix  whose  i-th  column  is  A(“. )  for 

x 


i=l,...,m.  Then  by  construction  Y =  {M1 x j x  6 { 0 , 1 } 
be  a  basis  for  Y.  For  i  =  we  find 


,  Let  Y(1>  ,...  ,Y{'^ 
such  that  M'x^  =  Y  ^ 


and  define  the  permutation  “j_=ITi 

by  construction  Y  = A ( < r, , . . . ,o „ > ) . 

1  i 


(i) 


( i ) 

-  m  .  Since  A(~.)  =  Y  ,  we  have 
m  l 


Now  let  M"  be  the  r*>.  boolean  matrix  whose  i-th  column  is 
A  (  -^ )  for  i=l,...,d..  Again  by  construction,  Y  =  {M"z|z  €  {0,1}  }. 
For  purposes  of  membership  testing,  it  suffices  to  have  the  list 
which  generate  the  coset  representatives  of  G/G  ^  . 


LEMMA  5.5.  The  set  RstOj1  ...  c  ^  |  z  , . . . ,  z.  €  { 0, 1} }  is  a  complete  set 
of  coset  representatives  for  G/G'1^. 


For  each  r  €  S  ,  let  ftJ(-)  =  (01  . . .  0  „  )  ~  where  M"z=A(~). 

n  k  i  x 

Clearly,  if  r  EG  then  f  (n )  €G^.  Hence  by  Lemma  2.1  we  have 

R 


LEMMA  5.6.  f  (RAND (G) )  «  RAND(G(1)). 

R 


Since  L =  (~  , . . . )  is  a  list  of  independently  chosen  elements  of  RAND(G), 
1  m 

o 

and  we  have  only  utilized  the  first  m  elements  of  L  to  construct  a  random 
presentation  of  G^  ^  ,  it  follows  from  Lemmas  2.2  and  5.6  that 


LEMMA  5.7.  L(1)  =  (fR(-mil),. 

chosen  elements  of  RAND(G^ 


The  above  Lemma  implies  that  we  can  repeat  the  above  construction,  to  construct 

(2)  .  ,  .  (1)  . 
a  random  presentation  of  G  from  the  first  m  elements  of  L  .A 

further  log  n  stages  yields  generators  for  the  entire  block-structure  tower 

G  =g('°)=>  G(1)  . .  .=>  G(h)  as  required.  The  linear  algebraic  computations 

(such  as  computing  basis  vectors) ,  required  in  each  stage  of  the  above 

ii  \  2.0(1) 

construction  of  G;  ,  can  be  done  in  time  0(log  n)  using  n  processors 

by  the  methods  of  [Borodin,  van  zur  Gothen,  and  Hopcroft,  82].  Since  there  are 
at  most  log  n  stages,  we  have: 

THEOREM  5.2.  Given  a  random  -presentation  of  2-group  we  cay:  construct 

aenerators  for  each  subgroup  of  the  block-structure  tower ,  in 
time  0(log  n)3  using  n° ^  processors. 


By  Theorems  5.1  and  5.2,  we  have: 


COROLLARY  5.1.  Given  a  random  presentation  of 
and  some  x  €  s  we  can  test  membership  in  G 


(worst  case)  2-group  G  cs^, 
in  expected  time  0(log  n)3 


processors. 
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COROLLARY  5.2.  Given  random  presentations  of  (worse  case)  2-group  G]_'G2-Sn 
we  can  test  in  expected  time  O(log  n)"3  using  processors. 


a 


T.  '--  Iv 
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